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Abstract: Coal mining has led to serious ecological damages in arid desert region of Northwest China. 
However, effects of climatic factor and mining activity on vegetation dynamics and plant diversity in this 
region remain unknown. Wuhai City located in the arid desert region of Northwest China is an industrial 
city and dominated by coal mining. Based on Landsat data and field investigation in Wuhai City, we 
analyzed the vegetation dynamics and the relationships with climate factors, coal mining activity and 
ecological restoration projects from 2000 to 2019. Results showed that vegetation in Wuhai City mostly 
consisted of desert plants, such as Caragana microphylla, Tetraena mongolica and Achnatherum splendens. And the 
vegetation fractional coverage (VFC) and greenness rate of change (GRC) showed that vegetation was 
slightly improved during the study period. Normalized difference vegetation index (NDVI) was positively 
correlated with annual mean precipitation, relative humidity and annual mean temperature, indicating that 
these climate factors might play important roles in the improved vegetation. Vegetation coverage and plant 
diversity around the coal mining area were reduced by coal mining, while the implementation of ecological 
restoration projects improved the vegetation coverage and plant diversity. Our results suggested that 
vegetation in the arid desert region was mainly affected by climate factors, and the implementation of 
ecological restoration projects could mitigate the impacts of coal mining on vegetation and ecological 
environment. 
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1 Introduction 


Vegetation is an important part of the global environment and affects several processes, including 
terrestrial carbon cycle, energy exchange and water balance (Cao and Woodward, 1998; Huang 
and Xue, 2016). The growth of vegetation is affected by climate factors, such as precipitation, 
temperature and wind speed. Vegetation is not only sensitive to climate but also regulates climate 
change through feedback of geophysical and biogeochemical mechanisms (Zhang et al., 2013). 
For example, Xu et al. (2018) analyzed the dynamic changes of vegetation in Mu Us Sandy Land 
from 1981 to 2013 using remote sensing technology and landform analysis methods, and found 
that the increase in vegetation was related to a reduction in wind power and vegetation 
colonization or afforestation on the bare sand dunes. The dynamic changes of Sahel vegetation 
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showed that increased precipitation was beneficial to vegetation restoration during 1984-2011 
using the Global Inventory Modeling and Mapping Studies (GIMMS-3g) dataset and field 
investigation methods (Dardel et al., 2014). Vegetation dynamics analyzed by the Moderate 
Resolution Imaging Spectroradiometer (MODIS) data of the Hulun Buir meadow steppe were 
mainly related to temperature, with no time lag, whereas the impact of precipitation was not 
obvious (Lin et al., 2015). In addition to climate factors, anthropogenic activities, such as open-pit 
coal mining, overgrazing and land use conversion from grassland to cropland, have destroyed the 
natural vegetation, leading to the degradation of ecological environment. For example, natural 
grassland was destroyed by cultivation, grazing activities and invasion of alien plant species (Wu 
et al., 2009). The mining has led to land desertification in Wuhai City of Inner Mongolia 
Autonomous Region, China (Ruan et al., 2020). In the Yushenfu coal mine in Shaanxi of China, 
the area of vegetation ecological degradation caused by coal mining was 42% higher than that of 
unexploited areas analyzed by normalized difference vegetation index (NDVI) data, indicating 
that coal mining activity had a negative impact on vegetation ecology (Liu et al., 2019). 
Eco-environmental projects could impact regional land use and vegetative coverage (Zhang et al., 
2018). Vegetation greenness has significantly increased over the past 27 years in eastern and 
southern America, owing to forest management interventions (Hicke et al., 2002). Vegetation area 
and NDVI in different climate zones of China from 1990 to 2015 increased by more than 50% and 
40%, respectively, after the implementation of ecological restoration projects (Shi et al., 2020). 

Vegetation fractional coverage (VFC) in Inner Mongolia showed an increasing trend from 1990 
to 2015. Spatially, VFC had an obvious decreasing trend from northeast to southwest in Inner 
Mongolia (Zhang, 2018). VFC in Inner Mongolia gradually decreases from east to west from 
1982 to 2010, and vegetation coverage was well correlated with precipitation (Tong et al., 2016). 
The area of reduced VFC was mainly distributed in the west of Alxa and scattered in the southern 
region in central Inner Mongolia from 2001 to 2010 (Tong et al., 2016). Wuhai City is an 
important coal production region in Inner Mongolia, China. Open-pit mining of coal in this region 
has caused a variety of environmental issues, such as land destruction, waste stacking, dust 
pollution, topsoil damage and soil degradation (Feng et al., 2019). Vegetation was degraded 
finally and ecological environment deteriorated (Dong et al., 2011). Deciduous shrub Tetraena 
mongolica Maxim that is a protected Zygophyllaceae plant distributed in Wuhai City (Wang et al., 
2014). Due to the contamination of heavy-metal molybdenum (Mo) and cobalt (Co) in the soil 
caused by coal mining, T. mongolica population was decreased (Wang et al., 2014). The spatial 
distribution of land damage in Wuhai City was featured by attenuation around the mining working 
face, a downward trend of NDVI from 2000 to 2009 and then an upward trend of NDVI from 
2009 to 2015 (Guo et al., 2018). With the continuous expansion of mining area and the increasing 
pressure of regional ecological environment, more attention needs to be paid to the vegetation 
dynamics in mining areas, typically in the arid regions. Over the past 20 years, the spatiotemporal 
changes in vegetation dynamics in the arid desert region of Northwest China were not well 
discussed, particularly in Wuhai City. The interactions between vegetation coverage and climate 
factors or human activities in the region remain unclear. Accordingly, the effect of coal mining on 
vegetation coverage and plant diversity in the regional scale was unknown. 

This study, based on Landsat TM (thematic mapper) and Landsat OLI (operational land imager) 
data, combined with field surveys, we analyzed vegetation dynamics and the relationships with 
coal mining activity and climate factors from 2000 to 2019 in Wuhai City. The aim of this study is 
to (1) explore the change in vegetation dynamics in Wuhai City from 2000 to 2019; (2) to analyze 
the relationship between vegetation and climate factors; and (3) to distinguish the effects of coal 
mining on vegetation coverage and plant diversity. This study could provide important 
information for ecological restoration in the coal mining area of arid desert regions. 


2 Materials and methods 


2.1 Study area 


Wuhai City is located in Inner Mongolia Autonomous Region, Northwest China (39°15’-39°52'N; 
106°36'-107°05'E; Fig. 1). The total area is approximately 1754.0 km?. Wuhai City is comprised 
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by three municipal districts, Haibowan district, Wuda district and Hainan district and the total 
population was 56.61x10° in 2019. The city belongs to the arid desert region with low 
precipitation, high evaporation and heavy sand-wind. The average annual precipitation ranges 
from 78.1 to 239.9 mm, with 90% occurring between May and September. The annual average 
temperature ranges from 9.2°C to 11.3°C, and the annual average relative humidity is 
approximately 37%—45%. The annual mean wind speed is approximately 2.3-3.1 m/s, and annual 
maximum wind speed is approximately 13.6—21.0 m/s. Sandstorms are frequent in this region, 
and strong winds mainly originate from northwest. Soils are mainly classified as gray desert soils 
and brown calcareous soils. The major vegetation types in the area include trees (Elaeagnus 
angustifolia, Pinus and Populus), shrubs (Artemisia ordosica, Haloxylon persicum, Caragana 
spp., Tetraena mongolica and Tamarix spp.), grasses (Stipa bungeana, Poa sphondyloides and 
Carex spp.) and other herbs (Agriophyllum) (Xu et al., 2018). 
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Fig. 1 Location of study area and distribution of sampling site in Wuhai City (a, b and c). Area M (d) is near to 
open-pit coal mine. Area M includes the core reserve area (M1), reserve area (M2) and marginal area (M3) of T. 
mongolica and open-pit coal mine (CM). BG, botanical garden; DR, dump waste reclamation area; TmD, T. 
mongolica disturbance area; W, wetland; CR, abandoned coking vegetation restoration area; SR, subsidence 
restoration area; NS, natural shrub vegetation area; TmN, T. mongolica Nature Reserve; GM, Gander Mountain. 


2.2 Field investigation 


Twenty-two sampling sites in Wuhai City were selected for field vegetation investigation in June 
2019. These sampling sites were divided into nine areas: T. mongolica Nature Reserve (TmN), T. 
mongolica disturbance area (TmD), natural shrub vegetation area (NS), Gander Mountain (GM), 
dump waste reclamation area (DR), abandoned coking vegetation restoration area (CR), 
subsidence restoration area (SR), botanical garden (BG) and wetland (W) near by the Yellow 
River. At each sampling site, three representative parallel quadrats were set up, and the latitude 
and longitude of the quadrats were recorded by a global positioning system (GPS). The size of 
each quadrat was 15 mx15 m. In addition, statistics were performed in each quadrat on vegetative 
indicators such as the coverage, height, plant density and composition. Alpha diversity indices 
such as Simpson's index, Shannon-Wiener's index and Pielou's evenness were calculated as 
described in Gonzdalez-Hernandez et al. (2020). Beta diversity represents the extent of species 
replacement or biotic change along an environmental gradient. Cody's index and Sgrensen's index 
of beta diversity were calculated as described in Gebrehiwot et al. (2019). 
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2.3 Remote sensing methodology 
2.3.1 Data sources 


In this study, Landsat 30-m resolution data were employed and obtained from the United States 
Geological Survey (http://earthexplorer.usgs.gov/), and a total of 13 stages of remote sensing 
images (Landsat TM or Landsat OLI) were selected from 2000 to 2019 (Table 1). To exclude the 
seasonal impact of crops within a year, this study focused on the growing season from May to 
August. In addition, images with few clouds and shadows covering the study area were 
preferentially selected when possible, or the influences were minimized. 


Table 1 Remote sensing images used in this study 


Remote sensing type Date (mm/dd/yy) Remote sensing type Date (mm/dd/yy) 
Landsat TM 08/29/2000 Landsat TM 06/18/2011 
Landsat TM 08/19/2002 Landsat OLI 07/28/2014 
Landsat TM 08/15/2003 Landsat OLI 05/30/2016 
Landsat TM 06/01/2005 Landsat OLI 05/17/2017 
Landsat TM 08/14/2006 Landsat OLI 08/29/2018 
Landsat TM 08/10/2007 Landsat OLI 06/19/2019 
Landsat TM 06/28/2009 
2.3.2 NDVI 


NDVI, calculated as the ratio of the difference between near-infrared and red visible reflectance 
values to their sum, is widely recognized as a method for studying the changes in terrestrial 
vegetation patterns and the ability of vegetation to absorb photosynthetically active radiation 
(Pettorelliand and Fang, 2005; Zhu et al., 2015). In this study, we calculated NDVI in ENVI 
(environment for visualizing images) v5.1 after removing the influence of clouds, atmosphere and 
solar elevation angle (Piao et al., 2001). Linear regression was applied to evaluate the strength of 
linear trends in NDVI over time. 

2.3.3 VEC 


VFC is a quantitative index for measuring vegetative cover, and we calculated VFC from NDVI 
image dataset using the dimidiate pixel model method. VFC is obtained as follows: 
VFC=(NDVI-NDVIsoi1)/(ND VIveg-ND VIsoit), (2) 
where NDVIveg is the NDVI of pure vegetation pixels; NDVIsoi is the NDVI of pure soil pixels; 
and VFC is the vegetation fractional coverage. Based on NDVI cumulative frequency table, we 
extracted a cumulative frequency of 5% for NDVIsoi and a cumulative frequency of 95% for 
NDV veg (Chu et al., 2019). The spatial distribution of VFC was then mapped across the study 
area. 
2.3.4 Greenness rate of change (GRC) 


GRC is defined as the slope of the smallest power-function linear regression equation for the 
inter-annual change in the synthetic NDVI during a given period. The least-squares method was 
used to simulate the temporal trends of average NDVI in the study area (Duan et al., 2011). The 
GRC is calculated as follows: 


nx >i x NDVI,- i} ix NDVI, 
GRC ee i=1 i=1 _i=1 , 
nx X P-Q iy 
i=1 i=1 
where GRC is the slope of this trend line; n is the number of years involved in the calculation; i is 
the position of the year involved in the calculation; and NDVI; is the average NDVI value of the 
first i year. The result reflects the trend of annual average NDVI during the study period. When 


GRC>0, NDVI values increase during the corresponding period; otherwise, NDVI values 
decrease (GRC<0) or remained constant (GRC=0). 


(3) 
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2.3.5 Land cover 

Land cover information is recognized as one of the crucial data components required for many 
aspects of global change studies and environmental applications (Sellers et al., 1995). For the land 
cover classification, raw remote sensing images were orthorectified and calibrated to 
top-of-atmosphere reflectance and used as the primary input data. We divided land cover scheme 
of Wuhai City according to the "Land Use Status Classification" national standard (GB/T 
2010-2017). Training sites with known classes were identified by field surveying, and their 
signatures were generated. Then, support vector machine (SVM) classifier was used for 
supervised classification of the image. SVM is a nonparametric supervised statistical 
classification technique. SVM-based techniques are quite practical for use in remote sensing 
because of their capability to produce higher classification accuracy even with limited training 
datasets (Mantero et al., 2005), and these techniques have been adopted for categorizing degraded 
coal mining land (Karan et al., 2016). Each pixel in the training dataset was further examined 
using high-resolution images in Google Earth™ and Landsat images to ensure that land cover 
labels were correctly assigned. In the accuracy assessment, a confusion matrix was created, which 
was represented by a table that shows correspondence between the classification result and a 
reference image assigned to a particular category. Then, the overall accuracy was calculated by 
the ratio of proportion of correctly classified pixels to total number of pixels in the confusion 
matrix. Kappa coefficient was also calculated to measure the agreement between classification 
and truth values. According to the evaluation, overall accuracy of each category classification 
exceeded 85%, and Kappa coefficient was above 0.7, indicating that final classification results 
had acceptable accuracy and small error. 


2.4 Climate factors 


Meteorological data were provided by the National Meteorological Bureau of China from 2000 to 
2019, including daily data of precipitation, temperature, relative humidity, mean wind speed and 
maximum wind speed observed from the meteorological station of Wuhai City. The annual 
precipitation (AP), annual mean temperature (MT), mean relative humidity (RH) and mean wind 
speed (WS) were calculated with the meteorological data. Strong wind frequency (W%) was 
calculated from the daily wind data and represented the percent of time during which the wind 
speed surpassed the threshold velocity required for sand entrainment. Thus we used 6 m/s at 10 m 
height as the threshold velocity in this study (Xu et al., 2018). 


2.5 Statistical analysis 


The Pearson's coefficients (r) were calculated between NDVI and climatic factors to investigate 
its possible relation. Ordinary least squares (OLS) regression was performed to detect the linear 
trends of NDVI and climate factors. Locally weighted scatter point smoothing (LOWESS) 
method was used to smooth the data of NDVI and climatic factors (Basarin et al., 2016). The 
slopes acquired from the OLS could quantify the changing rate, while LOWESS could depict the 
tendency in a finer timescale. 


3 Results 


3.1 Variations in VFC and GRC in Wuhai City 


VFC ranged from 0.00 to 1.00 in Wuhai City and was divided into five levels (Table 2). The 
temporal and spatial distributions of VFC showed obvious heterogeneity in Wuhai City from 2000 
to 2019 (Figs. 2 and 3). VFC in Wuhai City was mainly medium-high coverage (0.30—0.60) and 
high coverage (0.60—1.00) in 2000. From 2000 to 2009, low coverage (0.00—0.05) areas displayed 
an upward trend first and then a downward trend, while high coverage (0.60—1.00) areas 
continuously declined, indicating that vegetation degradation occurred during this period. 
Medium coverage (0.15-0.30) areas decreased dramatically, while high coverage (0.60—1.00) 
areas increased rapidly from 2009 to 2019, indicating that vegetation was improving at this stage. 
Moreover, the area with a high coverage (0.60-1.00) increased from 53.1% in 2000 to 60.9% in 
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2019, indicating that vegetation was improved during the 19 a period. These results suggested that 
vegetation degraded at first and then underwent a restoration process in the region. 


Table 2 Classification of vegetation fractional coverage (VFC) in Wuhai City from 2000 to 2019 
Level Value of VFC Description 


Moderate desertification land, rock, buildings, opencast mine, waste dump and 


Low coverage 0.00-0.05 low yield grassland 
Mediurn-low 0.05-0.15 Medium-low yield grassland, industrial land and sporadic vegetation 
coverage 
Medium coverage 0.15-0.30 Unutilized land, medium-low yield grassland and medium-low shrubland 
speed 0.30—0.60 Medium-high yield grassland, medium-high shrubland and sparse woodland 
High coverage 0.60-1.00 High-yield grassland, dense shrubland, reclaimed dump, cropland and dense 
woodland 
(a) 2000 (b) 2006 (c) 2009 (d) 2011 (e)2019 N 


VFC value 
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Fig. 2 Spatiotemporal distribution of vegetation fractional coverage (VFC) in Wuhai City. (a), 2000; (b), 2006; 
(c), 2009; (d), 2011; (e), 2019. 
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Fig. 3 Area changes in vegetation fractional coverage (VFC) in Wuhai City from 2000 to 2019 


We divided GRC into six categories based on the magnitude of the change (P<0.05; Table 3). 
Areas of moderate degradation (—0.10— —0.05) and slight degradation (—0.05-0.00) accounted for 
54.73% and 35.99% of total area in Wuhai City, respectively, indicating that vegetation degraded 
from 2000 to 2009. The area of moderate improvement (0.05—0.10) accounted for 83.10% of total 
area, indicating that vegetation improved from 2009 to 2019. The slight improvement area 
(0.00-0.05) accounted for 89.68% of total area from 2000 to 2019, and the area of vegetation 
improvement was larger than that of degradation, indicating that the overall vegetation generally 
improved during the 19 a period. 


3.2 Land cover in Wuhai City 


Land cover in Wuhai City was classified into ten categories (Fig. 4). The main water body area 
covered a branch of the Yellow River surrounded by wetland and cropland. Residential land was 
mainly distributed in the middle of Haibowan district. Coal mining area was mainly distributed in 
eastern and western part of Wuhai City, while industrial land was mainly distributed around the 
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coal mining area. Woodland was mainly distributed in the Gander Mountain. Grassland and 
shrubland were dispersedly distributed throughout the entire region. Barren land was sporadically 
found in the industrial area and around the coal mining area. 


Table 3 Dynamics of greenness rate of change (GRC) during three stages in Wuhai City 


2000-2009 2009-2019 2000-2019 
YS Category po ee Area MAT Atei P 
(xm?) of tota (km?) of tota (xm?) of tota 
area (%) area (%) area (%) 
< —0.10 Severe degradation 57.82 3.30 66.09 3.77 6.84 0.39 
—0.10- —0.05 Moderate degradation 960.02 54.73 24.08 1.37 22.56 1.29 
—0.05-—0.00 Slight degradation 631.30 35.99 26.01 1.48 102.29 5.83 
0.00-0.05 Slight improvement 49.26 2.81 121.68 6.94 1572.91 89.68 
0.05-0.10 „Moderate 17.97 1.02 1457.92 83.10 46.33 2.64 
improvement 
>0.10 Significant 35.86 2.04 58.23 3.32 0.18 0.01 
improvement 
(a) 2000 (b) 2006 (c) 2009 (d) 2011 (e)2019 y Legend 
Ka Barren land 
sa Grassland 
= Shrub land 
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= Woodland 
= Wetland 
= Water body 


= Coal mining area 
= Industrial land 
™ Residential land 


0 5 10km 
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Fig. 4 Dynamics of annual land cover classification in Wuhai City. (a), 2000; (b), 2006; (c), 2009; (d), 2011; (e), 
2019. 


The areas of water body, residential land, industrial land and coal mining area in Wuhai City 
significantly increased during the 19 a period, accounting for 3.7%, 11.4%, 2.9% and 10.6% in 
2019, respectively (Fig. 5). Water body area dramatically increased after completion of Wuhai 
Lake project in 2013 and encroached on the areas of barren land, grassland and shrubland 
thereafter. Increase in the areas of residential land, industrial land and coal mining were roughly 
in accordance with the rise in raw coal production after 2006 (Fig. 5). Areas of grassland and 
shrubland were occupied by the residential land, industrial land and coal mining, and vegetation 
in coal mining was degraded from shrubland to barren land. Areas of cropland and wetland were 
relatively stable, and might mainly be attributed to the 'cropland and wetland protection' policies 
of Chinese government. Areas of barren land tended to expand first and then decreased, while 
area of shrubland decreased first and then expanded moderately, and area of grassland markedly 
decreased after 2010. The afforestation and ecological restoration program were implemented in 
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Fig.5 Area of land cover and coal output in Wuhai City from 2000 to 2019 
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mining area, such as the backfilling of open-pit coal mine and reclamation of waste dumps, area 
of shrubland was improved after 2010. 


3.3 VEC in the coal mining area 


VFC in Wuhai City decreased first and then increased from 2000 to 2019, (Fig. 6a), the trend of 
VFC in area M (Fig. 1d) was in line with Wuhai City (Fig. 6a). VFC increased with the distance 
of the plots to open-pit coal mine (Fig. 6b). 

In 2000, VFC in area M revealed that medium-high coverage (0.30—0.60) and high coverage 
(0.60-1.00) was dominant, accounting for 20%-50% and 50%-80%, respectively (Fig. 6c). Area 
of low coverage (0.00-0.05) in M3 and CM plots gradually increased from 2000 to 2018, while 
area of high coverage (0.60—1.00) in M1 and M2 plots reached 98%, indicating that coal mining 
destroyed the land surface structure and dump waste occupied the T. mongolica area. Vegetation 
showed degradation in open-pit coal mine and surrounding areas during this period. The closer to 
the mine area a plot is, the lower the vegetation coverage and the more severe the impact of coal 
mining activity on vegetation is. 
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Fig. 6 Change in vegetation fractional coverage (VFC) over a chronological sequence in Wuhai City (a) and 
area M (b) and area percentage of different VFC levels in M1 (core reserve area), M2 (reserve area), M3 
(marginal area) and CM (open-pit coal mine) in 2000, 2009 and 2018 (c) 


3.4 NDVI in different sampling sites 


To reveal the vegetation dynamics in Wuhai City from 2000 to 2019, we analyzed NDVI at nine 
sampling sites (Fig. 7). NDVI in NS was increased over two decades, and provided evidence for 
the influence of climate factors on vegetation in the region. BG was built in 2003, vegetation in 
SR was restored with trees and shrubs from 2006, and the greening plants were irrigated by the 
drip irrigation system in both sites, and thus NDVI of the two sites was the highest. NDVI of GM 
was gradually increased under the influence of afforestation in recent years. TmN and W areas 
located in the nature reserves, and their NDVI was higher than that in NS. TmD was the closest to 
open-pit mine, and its NDVI was lower than that in TmN, indicating that vegetation was severely 
reduced by coal mining activity. NDVI of CR and DR was lower than that of TmN, 
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Fig. 7 Normalized difference vegetation index (NDVI) change at sampling sites in Wuhai City from 2000 to 
2019. BG, botanical garden; CR, abandoned coking plant restoration area; SR, subsidence restoration area; DR, 
dump waste reclamation area; TmD, T. mongolica disturbance area; NS, natural shrub vegetation area; GM, 
Gander Mountain; TmN, T. mongolica Nature Reserve; W, wetland. 
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suggesting that vegetation was destroyed by coal exploitation. Coal mining activity aggravated 
vegetation degradation, while the waste dumps reclamation and vegetation restoration in mining 
area had the potential to improve vegetation coverage. 


3.5 Plant community structure 


Plant community composition, coverage and plant diversity determined by field survey were 
shown in Table 4 and Figure 8. The value of VFC estimated by remote sensing analysis was 
highly correlated with the value of vegetation coverage measured by field survey (the decision 
coefficient reached 0.83; data not shown), indicating that VFC was a reliable indicator for 
evaluating changes of vegetation coverage at the regional scale. Vegetation in Wuhai City was 
mainly composed of arid desert plants, such as T. mongolica, Caragana microphylla, Salsola 
passerina, Zygophyllum xanthoxylon and Achnatherum splendens, according to field survey data 
from 2019 (Table 4). T. mongolica and Reaumuria soongarica were the dominant species in TmN. 
The number of T. mongolica decreased in TmD, while the number of Caragana microphylla and 
Zygophyllum xanthoxylon increased. Vegetation in GM was restored with trees in 2009 and 
dominated by Picea asperata and Pinus sylvestris, while Populus and Pinus sylvestris were 
dominant in coal mining SR area. The number of T. mongolica was the highest in TmN, followed 
by TmD, GM and NS2. However, T. mongolica was hardly found in BG, CR and DR (Fig. 8a), 
suggesting that native T. mongolica populations were damaged during coal exploitation and were 
difficult to restore under natural conditions. Vegetation coverage in TmN was higher than that in 
TmD, indicating that vegetation coverage was reduced by coal mining. DR had been restored by 
shrub species for 3 to 5 a, the vegetation coverage was higher than that of TmD. Vegetation 
restoration in CR was started in 2014 with trees and shrub species, vegetation coverage was 
similar to TmD. Vegetation greening by afforestation with trees in SR over 10 a, and thus its 
coverage was higher than those in other sites. Therefore, artificial restoration could increase 
vegetation coverage of the coal mining area. 


Table 4 Plant community composition in each sampling site 


Sampling site Species number Dominant species 
BG 15 Elaeagnus angustifolia; Populus alba var. pyramidalis; Platycladus orientalis 
SR1 14 Populus alba var. pyramidalis; Amorpha fruticose; Cynanchum chinense 
SR2 15 Pinus sylvestris var. mongolica; Picea asperata; Ulmus pumila 
CR 8 Elaeagnus angustifolia; Syringa vulgaris; Amorpha fruticosa 
DR1 9 Medicago sativa; Limonium bicolor; Amorpha fruticosa 
DR2 7 Artemisia salsoloides; Medicago sativa; Alhagi camelorum 
TmD1 12 Achnatherum splendens; Tetraena mongolica; Convolvulus tragacanthoides 
TmD2 7 Artemisia salsoloides; Convolvulus tragacanthoides; Achnatherum splendens 
TmD3 9 Caragana microphylla; Cynanchum hancockianum; Achnatherum splendens 
TmD4 10 Tetraena mongolica; Zygophyllum xanthoxylon; Achnatherum splendens 
TmD5 9 Tetraena mongolica; Zygophyllum xanthoxylon; Achnatherum splendens 
TmD6 7 Tetraena mongolica; Salsola passerina; Achnatherum splendens 
NS1 10 Achnatherum splendens; Ammopiptanthus mongolicus; Caragana microphylla 
NS2 11 Zygophyllum xanthoxylon; Nitraria sibirica; Cynanchum chinense 
GM1 11 Caragana microphylla; Picea asperata; Pinus sylvestris var. mongolica 
GM2 10 Caragana microphylla; Picea asperata; Pinus sylvestris var. mongolica 
TmN1 12 Tetraena mongolica; Achnatherum splendens; Reaumuria soongarica 
TmN2 11 Tetraena mongolica; Achnatherum splendens; Reaumuria soongarica 
W1 10 Phragmites australis; Sophora alopecuroides; Tamarix chinensis 
W2 10 Phragmites australis; Echinochloa crusgalli; Tamarix chinensis 


Note: BG, botanical garden; SR, subsidence restoration area; CR, abandoned coking plant restoration area; DR, dump waste reclamation 
area; TmD, T. mongolica disturbance area; NS, natural shrub vegetation area; GM, Gander Mountain; TmN, T. mongolica Nature 
Reserve; W, wetland. 
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Fig. 8 Plant alpha diversity indices, coverage, occurrence of T. mongolica (a) and beta diversity indices (b) at 
sampling sites in Wuhai City from 2000 to 2019. BG, botanical garden; CR, abandoned coking plant restoration 
area; SR, subsidence restoration area; DR, dump waste reclamation area; TmD, T. mongolica disturbance area; NS, 
natural shrub vegetation area; GM, Gander Mountain; TMN, T. mongolica Nature Reserve; W, wetland. 


Shannon-Wiener's index and Simpson's index among the sampling sites from low to high 
followed the order of TnD<NS<GM<W<TmN<BG indicating that species diversity in TmD was 
severely reduced by coal mining. The restoration age in CR and DR sites were shorter than that of 
SR, plant diversity and coverage in CR and DR were relatively lower than that in SR accordingly. 
The Pielou's indices in SR and TmN were nearly the same, but the dominant species were 
different. Cody's index increased from BG to DR (Fig. 8b), indicating that the high rate of species 
replacement occurred in coal mining restoration area. The highest variation in Sgrensen's index 
occurred between DR and NS sites, confirming the dissimilarity in plant communities between T. 
Mongolica area and restored areas (SR, DR and CR). Furthermore, the closer the coal mining area 
was, the lower the coverage and plant diversity was. These results suggested that ecological 
restoration projects might play important roles in accelerating ecological restoration and reducing 
the impact of coal mining on vegetation. 


3.6 Effects of climatic factors on NDVI 


From 2000 to 2019, NDVI in Wuhai City displayed a first decreasing and then increasing 
tendency (Fig. 9a). AP was about 199.1 mm and did not change significantly (Fig. 9b), while RH 
showed an increasing linear trend at a rate of 0.229%/a (Fig. 9c), indicating that the ecological 
environment was improved. MT first declined and then slightly increased, with an overall total 
negative linear trend at a rate of —-0.042°C/a (Fig. 9d). W% and WS declined significantly at rates 
of —0.018%/a and —0.032 m/(s-a), respectively (Fig. 9e and f), indicating that wind erosion risk of 
the study area had become less during the 19 a period. Overall, NDVI was positively related to 
RH (r=0.49, P=0.001), while NDVI was negatively correlated with W% (r= —0.67, P=0.012) and 
WS (r= —-0.32, P=0.007), suggesting RH and wind speed were important climate factors affecting 
vegetation activity in the region. 


4 Discussion 


4.1 Effects of climate factors on vegetation dynamics 


Wuhai City is a typical arid desert region, and vegetation is sensitive to climate change. 
Precipitation was the primary climatic factor that controlled vegetation growth in this region. 
NDVI gradually decreased from northeast to southwest in Inner Mongolia, showing an evident 
transition phenomenon between dry and wet areas (Li et al., 2016). NDVI was reduced to 
0.20—0.40 in the semi-arid steppe area in central area, indicating that vegetation in this area is 
gradually becoming sparser. In the arid desert region of western Inner Mongolia, NDVI is less 
than 0.10, with the surface mainly covered by desert, gobi, bare soil and bare rock (Li et al., 2016). 
Our result showed that AP in Wuhai City fluctuated around 199.1 mm during the 19 a period, the 


JOURNAL OF ARID LAND 


0.06 (a) Slope=0.0005, R?=0.379, P=0.025 144 (d) Slope= —0.0420, R?=0.151. P=0.090 


NDVI 


4001" 
99 (b) Slope= —0.1980, R?=0.001, P=0.924 
3004 
e : $ 
z s /\ 
E 200+ /\ 
a Aya fo = 4 e i 
= 1004 \7 e KPT 
v 
Fr T Ri y 
504 (c) Slope=0.2290, R?=0.419, P=0.002 
Z @ Jia- 2L 
9 Z a Sy a 
= nw 24 í 
n = 
14 
= = - - = T S T T T T T~ 
2000 2004 2008 2012 2016 2020 2000 2004 2008 2012 2016 2020 


Fig. 9 Inter-annual normalized difference vegetation index (NDVI, a), annual precipitation (AP, b), relative 
humidity (RH, c), mean temperature (MT, d), strong wind frequency (W%, e) and annual average wind speed 
(WS, f) in Wuhai City from 2000 to 2019. The linear trend (red dashed lines) is based on ordinary least squares 
regression, while the nonlinear trend (blue solid lines) is fitted by LOWESS (locally weighted scatter point 
smoothing). 


average NDVI of Wuhai City is less than 0.05, indicating that overall vegetation coverage was 
low, and the lower vegetation coverage was related to the lower AP. Precipitation most strongly 
and significantly limited vegetation growth over 45.1% of Inner Mongolia, while air temperature 
significantly limited vegetation growth over 0.73% at a small part of the northeastern Inner 
Mongolia (Tian et al., 2015). In Wuhai City, MT was increasing after 2012, and RH displayed an 
increasing trend though AP did not change during the 19 a period. In fact, the construction of 
Wuhai Lake was completed in 2013, and the area of water body reached 118.0 km?, drip irrigation 
facilities were established in plantation-style forests (GM) and shrubland in remediation sites 
(DR). Thus, the increased RH might be related to the extended water area and elevated soil water 
content caused by artificial irrigation. We infer that Yellow River Water Conservancy Program 
(2010) in Wuhai City may play an important role in improving air humidity. NDVI was positively 
related to RH, indicating that plant growth was mainly improved by soil water and RH in 2012. 
Soil losses caused by wind erosion decreased vegetation coverage and aggravated sand 
mobilization in the Sahel (Abdourhamane et al., 2019). Less wind could reduce the potential for 
surface erosion and sand transportation, and contributed to the growth and stability of vegetation 
in Mu Us Sandy Land (Xu, 2018). In Wuhai City, the reduction in wind speed and increase in RH 
over the 19 a period indicating that regional ecological environment was gradually improved. 


4.2 Effects of coal mining activity and ecological restoration project on vegetation 
dynamics 


Climate change was an important factor influencing vegetation in Wuhai City. However, 
anthropogenic activity was also an important driving factor. The dominant species of T. 
mongolica in Wuhai City significantly decreased with the strengthening of human disturbance 
from 1977 to 2005, such as mining, coking, roads, cutting firewood and over grazing (Zhen, 
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2012). T. mongolica population size and density were reduced dramatically by coal mining (Wang 
et al., 2014). From 2000 to 2009, VFC declined (Fig. 2), with GRC<0 (Table. 3) and vegetation 
was seriously degraded, especially in the coal mining area (Fig. 7). One potential explanation was 
that vegetation growth was inhibited by the low precipitation and temperature in this period (Fig. 
9), and the other was that vegetation was destroyed by increasing industrial land and coal mining 
activity (Figs. 4 and 5). From 2009 to 2019, VFC in Wuhai City increased, with GRC>0, and 
vegetation showed an upward trend, particularly in the coal mining area (Fig. 7). Further, the 
trend of vegetation coverage showed diversity in local and regional scales, and different land use 
types had different changes in vegetation coverage (Fig. 7), indicating that meteorological factors 
and anthropogenic activities might act together to affect the spatiotemporal changes of vegetation 
coverage. Our field investigation found that number of T. mongolica, coverage and plant diversity 
in the coal mining area was decreased compared with that in T. mongolica Nature Reserve. The 
closer to the coal mining area the site is, the lower the plant diversity is. This could be caused by 
land occupation and heavy metal contamination in the soil (Wang et al., 2014). 

Facing vegetation degradation and land desertification, a series of national ecological 
restoration projects or programs were implemented to restore degraded ecosystems in Wuhai City, 
including Beijing—Tianjin Wind/Sand Source Control Program (2000), the fourth phase of 
Three-North Shelterbelt Development Project (2001), Grain for Green project (2001), Wildlife 
Protection and Nature Reserve Development Program (2001) and Grazing Ban (2003). Based on 
these projects, government established West Erdos Natural Reserve to conserve T. mongolica with 
an area of 139.1 km? in Wuhai City. NDVI in TmN area gradually increased (Fig. 7), and area of 
shrubland in Wuhai City increased after 2008 (Fig. 5). Afforestation practice was initiated in 
barren mountains and SR1 of Wuhai City in 2006, with greening area covering over 3.3 km’. BG 
was constructed in 2003, covered an area of 1.4 km? with 1.1 km? of green land. Meanwhile, 
construction of Gandel Mountain Ecological Civilization Scenic Spot was initiated, and 6.7 km? 
of mountain afforestation (GM) was completed in 2009. Area of barren land had been reduced 
since 2006 (Fig. 5), and NDVI in study area had been increased since 2009 (Fig 2). Thus, 
greening increased vegetation coverage and could further reverse land desertification. Evident 
land greening phenomena in sand and desert areas of the Loess Plateau of China (including Wuhai 
City) were observed during 2006-2008, and NDVI rapidly increased during 2011-2015 (Liu et al., 
2020). NDVI also significantly increased in Inner Mongolia from 2000 to 2012 (Tian et al., 2015). 
In addition, China's greening was attributed to large-scale tree plantations in low-productive areas 
from 2000 to 2017 (Chen et al., 2019). These previous studies coincided with our results. 
Moreover, Wuhai City also implemented several ecological construction policies and projects in 
order to improve the ecological environment in the coal mining area, such as Opinions on 
Strengthening the Supervision and Administration of Mineral Resource Development project 
(2008), Ecological Construction project (2012) and Comprehensive Environmental Treatment 
project (2014). Accordingly, revegetation in waste dump was restored with native shrub and herb 
species in 2014. Roads in the coal mining area and abandoned cocking plant site were afforested 
with tree and shrub pattern. Consequently, NDVI in the coal mining area showed a significant 
improvement in recent years (Fig 7). Revegetation of waste dump also occurred in open-pit coal 
mining area in the past decades, and its vegetation coverage has reached more than 90% (Wang et 
al., 2020). Therefore, vegetation restoration also greatly improved ecological environment in the 
regional scale. Overall, NDVI in Wuhai City showed a consistent increasing trend in recent years, 
while potential degradation risks still existed in parts of the coal mining area. Anthropogenic 
impacts that are superimposed on the long-term climate-driven variations are also significant for 
vegetation restoration at both local and regional scales (Xu et al., 2018). However, the 
contribution of human activities to vegetation greening was not yet quantitatively assessed in this 
region. Therefore, plant diversity and community structure in the restoration area need to be 
further examined. Moreover, ecological restoration projects should fully consider regional 
resource endowments and carrying capacity of the environment. 
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5 Conclusions 


We revealed vegetation dynamics of Wuhai City in an arid desert region based on the land 
satellite remote sensing data and field survey. Vegetation in Wuhai City was mainly consisted of 
desert plant species, such as Caragana microphylla, T. mongolica and Achnatherum splendens. 
From 2000 to 2009, vegetation degradation in Wuhai City was induced by both climate factors 
and coal mining activity, climate factors might be the main factors, and coal mining accelerated 
the vegetation degradation. From 2009 to 2019, the increase in vegetation was attributed to the 
implementation of ecological restoration projects and improvement of ecological environment. 
Overall, climate was the main reason for vegetation dynamics, coal mining affected vegetation 
coverage and plant diversity, while ecological restoration policies and projects played important 
roles in vegetation restoration. The study provided reference for ecological restoration in the arid 
desert regions. 
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